The Small Polaron Crossover: Role of Dimensionality 
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The crossover from quasi free electron to small polaron in the Holstein model for a single electron is 
studied by means of both exact and self-consistent calculations in one dimension and on an infinite 
coordination lattice, in order to understand the role of dimensionality in such a crossover. We 
show that a small polaron ground-state occurs when both strong coupling (A > 1) and multiphonon 
(a 2 > 1) conditions are fulfilled leading to different relevant coupling constants (A) in adiabatic and 
(a 2 ) anti adiabatic region of the parameters space. We also show that the self-consistent calculations 
obtained by including the first electron-phonon vertex correction give accurate results in a sizeable 
region of the phase diagram well separated from the polaronic crossover. 



I. INTRODUCTION 

Recent optical measurements of the insulating par- 
ent compounds of the high-temperature superconductors 
show the presence of polaronic carriers , and evidence 
for intermediate and strong lattice distortions has been 
given also for the colossal magnetoresistance manganitcs 
Q and Nickel compounds [||. The recent observation of 
one-dimensional stripes in the high-temperature super- 
conductors Q and in manganites suggests a comprehen- 
sive study of the role of dimensionality in the polaronic 
crossover. A detailed study of the small polaron crossover 
is demanded also by the recent experimental results on 
manganites ||. 

The polaronic state is characterized by strong local 
electron-lattice correlation and is a non-perturbativc phe- 
nomenon. It therefore cannot be described by simple 
summation of the perturbative series such as the one 
which defines the Migdal-Eliashbcrg (ME) theory. Here, 
we provide a detailed study of the crossover which occurs 
at intermediate electron-lattice couplings from quasi-free 
electron to small polaron ground state, with a particular 
emphasis on the role of system dimensionality. 

We consider the simple Holstein molecular-crystal 
Hamiltonian for a single electron, which reads: 

H = -t 4°] + .9 4°i (°>i + 4) + w o Y 4 a% w 

(ij) i • 

where c, (ej) and oti (oJ) are, respectively, the destruction 
(creation) operators for an electron and for a dispersion- 
less phonon of frequency uiq on site i. The hamiltonian 
(Q) represents a non-trivial many-body problem and it 
has been already studied in recent years by means of nu- 
merical and analytical JTo| - |l2[ techniques. 



Two dimensionless parameters are introduced to mea- 
sure the strenght of electron-phonon (el-ph) interaction: 
A = g 2 /(Dujo) and a — g/ujo, where D = 2td is the 
half-bandwidth for the free electron and d is the system 
dimensionality. 

A is originally introduced in the weak coupling pertu- 
bation theory (g/t <C 1) and is the coupling parameter 
of a ME approach in the case of one electron. It can also 
be viewed as the ratio between the small polaron energy 
Ep = —g 2 /ujQ and the free electron energy -Ef ree = —D. 

The parameter a is the relevant coupling in the atomic 
limit (t = 0). In this limit a measures the lattice dis- 
placement associated to the polaron and a 2 is the aver- 
age number of phonons bound to the electron. According 
to the Lang-Firsov results followed by the Holstein ap- 
proximation, a also rules the reduction of the effective 
hopping t* = iexp (-a 2 ) f§|p|. 

Besides A and a, the el-ph system described by eq.(0) is 
governed also by another dimensionless parameter: LOo/t. 
It measures the degree of adiabaticy of the lattice mo- 
tion (lattice kinetic energy ~ u>o) compared to the elec- 
tron one (electron kinetic energy ~ t). In the adiabatic 
regime (u>o/t C 1), A > 1 is a condition sufficient to 
give a polaronic state since the electron is bound to the 
slowly moving lattice giving rise to a strong enhancement 
of effective mass. In the antiadiabatic regime (u>o/t ^> 1) 
such a picture is no longer true due to the fast lattice mo- 
tion. In this case, polaronic features such as strong local 
electron-lattice correlations arise only when the electron 
is bound to a large number of phonons (a 2 > 1). To 
summarize, in both adiabatic and antiadiabatic regimes, 
a polaronic state is formed when both A > I and a 2 > I 
inequalities are fulfilled ||] . This conclusion is in contrast 
with ref. |l(J where it is argued that A > 1 is the only 
condition for small polaron formation. 
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FIG. 1. Ground state energy results in d=l. The exact 
diagonalization results are compared with the NCA (short 
dashed) and VCA (long dashed) calculations. The arrows 
mark A c (see Fig.||). 
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FIG. 2. Ground state energy results for an infinite coor- 
dination lattice. Comparison between dynamical mean field 
(solid line), NCA (short dashed) and VCA (long dashed). The 
arrows mark A c (see Fig.H). 



The parameter ivg/t influences also the dependence of 
the behavior of the el-ph coupled system on the sys- 
tem dimensionality. We shall show that in the antia- 
diabatic regime the small polaron formation does not de- 
pend on the system dimensionality On the other hand, 
dimensionality plays a crucial role in the adiabatic regime 
uiq <C t. This can be traced back to the adiabatic limit 
uJo/t = 0. In fact, in d=l the ground state is localized 
for any finite value of A and a crossover occurs between 
large and small polaron at A ~ 1, whereas for d > 2 it has 
been shown that a localization transition occurs at finite 
A from free electron to small polaron |Q. The different 
adiabatic behaviors between Id and 2d systems could be 
relevant to describe the motion of polarons as defects on 
top of Id charge striped structures such as those observed 
in cuprates 0] and manganites Pi- 



ll. RESULTS 

We study the relevance of u>o/t and of the lattice 
dimensionality d by using two alternative exact cal- 
culations: exact diagonalization of small one dimen- 
sional clusters (ED- Id) and dynamical mean field the- 
ory (DMFT-3d). In the ED-ld approach, the infinite 
phonon Hilbert space must be truncated to allow for a 
given maximum number of phonons per site n m ax- In 
order to properly describe the multiphonon regime (ex- 
pecially in the adiabatic regime where a large number 
of low energy phonons can be excited) we chose a cut- 
off of Umax = 20. This high value forced us to restrict 
our analysis to a four-site cluster in the strong-coupling 
adiabatic regime. 



In the weak-coupling regime and for larger phonon fre- 
quencies a lower value of nmax is needed, allowing us 
to consider larger clusters up to ten-twelve sites. We 
checked that finite-size effects do not affect the crossover 
coupling, since small-polaron formation is a local, high 
energy process. 

The dynamical mean field theory approach can be seen 
as the exact solution of the small polaron problem on an 
infinite coordination lattice |ll[| . The formulation of the 
DMFT requires the knowledge of the free particle DOS. 
A semi-circular DOS can mimic a three-dimensional case: 
in the following we will therefore refer to this approach 
to as DMFT-3d. 

We calculate the exact ground state energy Eq ob- 
tained by means of ED-ld and DMFT-3d and we com- 
pare the results with the self-consistent non-crossing 
(NCA) and vertex corrected approximations (VCA). 
These two approximations are defined by the self- 
consistent calculation of the electronic zero-temperature 
self-energy to) given below: 

E(fc > u) = 2^ a *EpG(p J w-wo) x 
l + ^E q G(q-p + k,LO-u Q )G(q,LU-2LU )~\ , (2) 

where G(k,uj) is the retarded fully renormalized single 
electron Green's function: 



G(fc, (jj) 1 = UJ — £fe — S(fc, u) + iS. 



(3) 



which will be determined self-consistently. From 
eqs.(^||), the ground state energy E is given by the low- 
est energy solution of Re G(k, E Q )~ l = 0. The NCA ap- 
proach amounts to compute S by retaining only the 1 in 
the square brackets of eq. (g) . The VCA is instead given 
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by the inclusion also of the second term in square brack- 
ets of eq.(H) which represents the first vertex correction. 
This approach is formally similar to the approximation 
scheme used in the formulation of the non-adiabatic the- 
ory of superconductivity [ p"5[ and a comparison with ex- 
act results therefore provides also a test of reliability of 
such an approach for the one-electron case. 

In figjl] we compare the ground-state energy Eq ob- 
tained by ED-ld with the NCA and VCA results. The 
same quantities evaluated in the DMFT-3d case are 
shown in fig. [2| We have chosen the same half-bandwidth 
D in both DMFT-3d and ED-ld. In the adiabatic regime, 
the agreement of both approximations with exact results 
strongly depends on the system dimensionality as a re- 
sult of the different low-energy behaviour of the DOS . 
In fact, moving from (Jo/t = 0.2 to (Jo/t = 0.5, before 
the crossover the agreement of the self-consistent calcu- 
lations with the exact results is improved for the Id case 
(fig. |l|) whereas it becomes poorer for the 3d one (fig. 
||). However, the VCA approach represents a significative 
improvement with respect to the NCA for every system 
dimensionality and over a significant range of parameters. 

As it is seen from the comparison of fig. |l| and |[ for 
large u>o/t both approximate and exact results tend to 
become independent of dimensionality. This can be un- 
derstood by realizing that in this regime the system can 
be thought as a flat band "atomic" system in interaction 
with high energy phonons. It is also clear from figs. |l| 
and ^ that both the self-consistent NCA and VCA calcu- 
lations deviate from the exact results when the crossover 
towards the small polaron regime is approached. 

A complete comparison between the exact results and 
the VCA approach in the parameter space X-LOo/t is 
shown in figs. |^(a)-(b). We explicitly evaluated both 
in Id (a) and 3d (b) the relative difference 5Eq — 
2 |£yCA_£exact| / |£yCA + £Cxact| where ^exact and 

E^^^ are the ground-state energies evaluated by exact 
techniques and the VCA, respectively. To analyze the re- 
gion in the parameter space where the VCA agrees within 
a given accuracy with the exact results we report in figs. 
||(a)-(b) lines of constant SEq. 

As already mentioned, the agreement between self- 
consistent approximations and exact results is sensible 
to dimensionality. For d > 2, approaching the adiabatic 
limit and for small couplings the electron tends to be 
free. For this reason self-consistent approximations work 
well. On the contrary, in the adiabatic limit and for d=l 
the ground state is a localized large polaron and self- 
consistent approximations fail to predict its energy. In 
general, VCA (and so NCA) works well outside the po- 
laron region whatever polarons are either small or large. 
This can be seen directly from figs. ||(a)-(b) where the 
critical coupling Ac of the crossover to small polaron is 
depicted as a dotted line. The critical coupling A c is de- 
fined as the value at which dE^jdg has maximum slope. 



In the same figures, we provide also an estimate of the 
width of the crossover (shaded areas) obtained by looking 
at the maximum slope of \d 2 Eo/dg 2 \. We checked that 
different criteria, like e.g. the effective mass enhancement 
pT|, provide the same qualitative results. 
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3. Phase diagram in the X-uo/t plane 



one- dimensional (a) and the infinite coordination lattice (b) 
Holstein model. The dotted line is the polaron crossover value 
A c and the width of the crossover is evidentiated by a shaded 
area. The isolines represents the relative difference between 
the exact and the VCA result for the ground state energy (see 
text) . 

III. CONCLUSIONS 

In conclusion, we have shown that the crossover to- 
wards the small polaron state depends strongly on the 
adiabaticity parameter uJo/t. In the antiadiabatic regime 
the crossover is ruled by a 2 and it is independent of the 
system dimensionality. In the adiabatic regime the rel- 
evant coupling is A and the crossover occurs from large 
to small polaron in Id, while in 3d the crossover is from 
quasi free electrons to small polarons. In the latter case 
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self-consistent approximations work better than in Id 
systems. We have also shown that self-consistent cal- 
culations provide ground state energies which agree well 
with exact results outside the small and large polaron re- 
gion of the phase diagram and that such an agreement is 
increased when vertex corrections are taken into account. 
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